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Abstract. The steady state packing fraction of a tapped granular bed is studied for 
different grain shapes via a discrete element method. Grains are monosized regular 
polygons, from triangles to icosagons. Comparisons with disk packings show that the 
steady state packing fraction as a function of the tapping intensity presents the same 
general trends in polygon packings. However, better packing fractions are obtained, as 
expected, for shapes that can tessellate the plane (triangles, squares and hexagons). In 
addition, we find a sharp transition for packings of polygons with more than 13 vertices 
signaled by a discontinuity in the packing fraction at a particular tapping intensity. 
Density fluctuations for most shapes are consistent with recent experimental findings 
in disk packing; however, a peculiar behavior is found for triangles and squares. 
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1. Introduction 

Granular materials settle under gravity and come to mechanical equilibrium unless an 
external excitation is provided. The properties of such static packings are difficult to 
predict, since the history of preparation of the sample is important. However, there exist 
different protocols to prepare a granular bed in a well defined macroscopic state. In such 
state, the packing fraction (and other macroscopic observables such as the pressure on 
the container) are reproducible if the given protocol is followed. A canonical example 
of this is the steady states obtained by tapping the sample with a given intensity [JJ. 
After a suitable annealing, tapping at a constant intensity produces mechanically stable 
configurations (inherent states, or microstates) whose ensemble has well defined mean 
values of all macroscopic observables. 

In recent years, the dependency of the steady state packing fraction, 0, on the 
tapping intensity, T, has been shown to be nonmonotonic; presenting a minimum at 
relatively high values of T for disks and spheres [21 G], and a maximum at very low V for 
spheres [TJ. In general, the symbol T is used for the reduced peak acceleration given to 
the system during a tap. However, we will use T in what follows to refer to any suitable 
parameter that characterizes the tapping intensity. 

On the one hand, there exist some studies on the response to tapping of non- 
spherical particles El [71 IE] , however these do not consider polygonal particles. On 
the other hand, the are some investigations on polygon packings [TO], [11]. These 
latter studies, however, do not focus on the steady state obtained after a repeated pulse 
excitation. Inspired by previous works on pentagon packings [ToT 1TB] . we investigate 
the (fy-T tapping curve in the steady state for monosized regular polygons with different 
number N of vertices; from triangles (N = 3) to icosagons (N = 20). As the number of 
vertices grows, we expect polygon packings to approach the properties of disk packings. 
Since depending on the number of vertices these particles may or may not tessellate the 
plane, we also expect strong deviations from the general trends for some grain shapes. 

In this paper, we compare the general features found in the <p-T curve of disk 
packings with those of regular polygons. Although some general trends are conserved, 
new phenomenology emerges. 

In Section [2] we present the simulation technique and the model particles. In Section 
I3.1l we analyze the behavior of polygons with fewer than ten vertices. In Section [372] we 
present results for polygons of up to twenty vertices. Section 13.31 is devoted to the study 
of the density fluctuations. Finally, we draw the conclusions in Section H] and point out 
some interesting areas of research suggested by the new results. 

2. Simulation 

We perform molecular dynamic type simulations by solving the Newton-Euler equations 
of motion for rigid bodies confined on a vertical plane. Gravity acts on the negative 
vertical direction. The bodies (particles) are placed in a rectangular box which is 
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confined to move in the vertical direction. This box is high enough to avoid particles 
to contact the ceiling during the simulations. We prepare nineteen samples that consist 
of 500 monosized regular polygons of a single type (from triangles to icosagons) or 
monosized disks. Particles, initially placed at random without overlaps in the box, are 
let to settle until they come to rest in order to prepare the initial packing. Then, the 
same tapping protocol is applied to each sample. 

We set the particle-particle interactions to yield a normal restitution coefficient 
e = 0.058 and a static and dynamic friction coefficient fi s = fid = 0.5. The confining 
box is 24. 8r wide and 2000r tall (with r the radius of the particles). The particle- 
box friction coefficient is fi s = fid = 0.07 and the restitution coefficient is as in the 
particle-particle interaction. All polygons have the same radius and material density. 
Therefore, the actual weight of a particle depends on the number of vertices. We use as 
unit mass, m, the mass of a disk; as unit length r; and the unit time is (r/g) 1 ^ 2 , with g 
the acceleration of gravity. 

Tapping is simulated by giving the box an impulse. In practice, we set the initial 
velocity v of the box (originally at rest after deposition) to a given positive value and 
restart the dynamics. In doing so, the box and its filling move upward and fall back on 
top of a zero restitution base. While the box dissipate all its kinetic energy on contacting 
the base, particles inside the box bounce against the box walls and floor until they fully 
settle. After all particles come to rest a new tap is applied. The intensity of the taps is 
measured by the initial velocity imposed to the confining box at each tap (i.e. T = vq). A 
similar parameter (the lift-off velocity) has been recently proposed as a suitable measure 
of the tap intensity [T7] . 

The tapping protocol consist in a series of 50000 taps. Every 250 taps we change 
the value of T by a small amount AI\ We initially decrease T from « i5.0(rg) 1,/2 down 
to a very low value and then increase it back to its initial high value. At each value of 
T the last 150 taps are used to average the packing fraction in order to plot the cfy-T 
curve. 

The simulations were implemented by means of the Box2D library [18J. Box2D uses 
a constraint solver to handle hard bodies. At each time step of the dynamics a series 
of iterations (typically 20) are used to resolve penetrations between bodies through a 
Lagrange multiplier scheme [TH]. After resolving penetrations, the inelastic collision at 
each contact (a contact is defined by a manifold in the case of polygons) is solved and 
new linear and angular velocities are assigned. The equations of motion are integrated 
through a symplectic Euler algorithm. The time step 5t used to integrate the equations 
of motion is 0.025-Jd/g. Solid friction is also handled by means of a Lagrange multiplier 
scheme that implements the Coulomb criterion. This library achieves a high performance 
when handling complex bodies such as polygons. 



Steady state of tapped granular polygons. 



4 




Figure 1. Mean packing fraction <f> as a function of tapping intensity T for triangles 
(violet), squares (red), pentagons (green), hexagons (blue), heptagons (yellow), 
octagons (cyan), nonagons (magenta), and disks (black). Except for disks, all curves 
correspond to a progressive decrease of T followed by an increase back to high values. 
For disks only the decreasing part has been carried out. Error bars correspond to the 
estimated error of the mean. 

3. Results 

3.1. From triangles to nonagons 

The steady state packing fraction as a function of the tapping intensity for triangles, 
squares, pentagons, hexagons, heptagons, octagons and nonagons is presented in Fig. [T] 
alongside with the results for disks. Packing fraction is estimated from the number 
density measured in a rectangular slab of half the packing hight at the middle of the 
sample. The fact that the same (p-T curve is obtained for decreasing and increasing T 
indicates that these states are reversible and that is uniquely defined for each T. From 
this results we can see that polygon packings present similar features to those observed 
in disk packings. At low tapping intensities, a decrease of <p is observed for increasing T 
down to a minimum packing fraction (j) m i n . A further increase of T induces an increase 
of (f) until a plateau is reached at a packing fraction somewhat lower than the maximum 
obtained for the lowest values of V . Disks also show a not very pronounced maximum 
at low T which is not observed in polygon packings. This maximum has been recently 
observed in sphere packings [3]. Another overall trend is that the range of packing 
fractions attained by disk packings is narrower than for polygons. 
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Figure 2. Minimum packing fraction 4> m in as a function of the number of vertices. 
The blue line is drawn only to guide the eye. 



Beyond these general features, there are some peculiarities associated to the ability 
of a given polygonal shape to tessellate the plane. As is to be expected, triangles, squares 
and hexagons can reach packing fractions of nearly 1 at the lower tapping intensities. All 
other shapes reach packing fractions similar to disk packings at low V. It is important to 
notice at this point that our results differ from those obtained by Vidales et al. [IHJ [16] 
in the case of pentagons in the framework of a pseudo-dynamic algorithm. In Refs. 
[T5| 116] the (jy-T curve does not present any minimum of the packing fraction. 

In Fig. [2j we plot the minimum steady state density, (fi m in, as a function of 
the number of vertices of the polygon. As the number of vertices is increased, a 
consistent increase of 4> m i n is found for all polygons with the exception of triangles, 
squares and hexagons. As we mentioned, these three polygons can tessellate the plane. 
Correspondingly, triangles, squares and hexagons present higher densities than expected 
by the trend showed by all other polygons. We have seen that the position, T min , of 
the minimum is independent of the number of vertices. The existence of <f) m i n has been 
associated to a competition between arch formation and arch breaking [2] . The position 
T min of such minimum signals the crossover between a regime where arches cannot form 
due to the particles settling one by one (in a sequential manner) at very high T, and 
a regime where arches do form but are "melted down" in successive taps creating a 
dynamic equilibrium. The fact that T m i n is the same for all shapes is a clear indication 
that arching is not favored (nor prevented) by any particular shape at these intermediate 
values of T. 
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Figure 3. Mean packing fraction ^ as a function of tapping intensity T for nonagons, 
decagons, ... and icosagons. The brown data in the lower right panel correspond to 
disks. A progressive decrease (red data) of T is followed by an increase (blue data) 
back to the high initial values. Error bars as in Fig. [T] 



3.2. An unforeseen sharp transition for triskaidecagons and beyond 

We now focus on the behavior of polygons with larger number of vertices (from nonagons 
up to icosagons). Figure [3] shows the (fi-T curves for each shape. One might have 
expected that a smooth change would appear in these curves as the number of vertices 
is increased up to a point where the behavior of the n-vertex polygon will converge to 
the one shown by disk packings. However, a sudden change is found as we move from 
dodecagons to triskaidecagons. While a continuous (f>-T curve is observed for polygons 
with up to 12 vertices, a sharp discontinuity in <fi is present in all packings with polygons 
of 13 vertices or more. A gap of "forbidden" values of 4> appears between roughly 0.80 
and 0.83 in all these polygon packings with more than 12 vertices. It is important 
to mention that fluctuations are rather large, and configurations (microstates) with 
0.80 < < 0.83 are rather common. It is the mean values that present a gap. 

A similar discontinuity has been seen in tapped disk packings simulated under 
a pseudo-dynamic algorithm [20J. However, this is not observed in our simulations 
of disks (see brown data in the lower right panel in Fig. [3]) nor in previous molecular 
dynamic simulations where the same region of T was explored [21] . The pseudo-dynamic 
algorithm [20] conducts a deposition of disks that roll on top of each other without 
sliding. This might mimic, rather realistically, the behavior of regular polygons with 
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Figure 4. Mean packing fraction as a function of tapping intensity T for 
tctrakaidccagons. Panel (a), increasing T. Panel (b), decreasing V. The red and 
blue data correspond to independent realizations of the tapping protocol. The black 
data correspond to the ones presented in Fig[3]for tetrakaidecagons, where larger steps 
in T arc taken. The full and dashed black lines are to guide the eye. Error bars as in 
Fig.ffl 



a large number of vertices. These polygons behave like gears in the sense that they 
interlock very easily just as if they were infinitely rough disks. We presume this basic 
characteristic shared by polygons with many vertices and disk that roll without sliding is 
the underlying phenomenon that leads to the emergence of a discontinuous <p-T curve. 
We mention in pass that, although it is difficult to relate with the static packings 
studied here, a similar discontinuity has been reported in an oscillation experiment of a 
2D granular sample [22] . 

In order to have a rough indication of the nature of the transition, we have made 
a more detailed simulation for tetrakaidecagons (N = 14). In Fig. HI the steady state 
value of is plotted for T in the interval [2.8,4.0] with a smaller Ar step. In panel 
(a), we plot two independent experiments obtained by increasing T alongside with the 
corresponding results from Fig. [3] (where a larger Ar was used). The results for the 
reversed protocol in which F is decreased is presented in panel (b) of Fig. HJ In Fig. Ufa), 
the system seems to present a first order type transition where metastable branches are 
explored. Since fluctuations are rather large for this small system sizes, the system 
may explore microstates compatible with both "coexisting" phases. Nevertheless, in 
Fig. Mb), where the protocol corresponds to decreasing T, the transition looks much 
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Figure 5. Standard deviation A(j> of the packing fraction in the steady state as a 
function of <fi. The red line is a simple running average to guide the eye. The arrows 
indicate the direction of increasing T. 



smoother if the rate Ar is reduced. Although the data is noisy, we can see that the 
width of transition region is rate dependent. 



3.3. Density fluctuations 

Density fluctuations have recently received renewed interest as a way to measure 
configurational temperature (as defined by Edwards [23]) and entropy [21]. It was 
in a fluidization experiment that a nonmonotonic dependence of the fluctuations A0 as 
a function of in the steady state was first reported [22] • in that work, Schroter et al. 
found a minimum in the density fluctuations for spheres. However, a recent study on 
disks reported a maximum in fluctuations from both, experiments and simulations [3]. 

In Fig. [2] we show the steady state density fluctuations A0 as measured by the 
standard deviation as a function of for several polygons and disks. The results for 
disks are entirely in agreement with Ref. [3]. A clear maximum in A0 appears for 
disks. One can also see that states of equal at each side of (fi m in present slightly 
different fluctuations. This indicates that these states are not equivalent and that is 
not sufficient to characterize the macroscopic state. A more detailed analysis of this 
can be found in Ref. [3] where the force moment tensor is found to be a suitable extra 
macroscopic variable in accordance with theoretical suggestions [26J. 

The behavior of the density fluctuations in the polygon packings show the signal of 
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the transition for shapes with iV > 12. However the same general trends as those seen 
for disks are observed. Interestingly, a peculiar behavior appears for triangles, squares 
and pentagons. Pentagons present the same fluctuations at both sides of the minimum 
whereas triangles and squares present a reversed situation where fluctuations are larger 
for large T, instead of smaller as seen in all other shapes. This change in trend should 
have an important impact in the calculation of configurational temperature and entropy. 
We will pursue this point further elsewhere. 

4. Conclusions 

We have carried out simulations of the tapping of assemblies of regular polygonal grains 
and studied the steady state of such systems. The comparison with more widely 
studied disk packings has shown some general similarities but also remarkable new 
phenomenology. 

On the one hand, beyond the expectable result for triangles, squares and hexagons 
that cover the space if gently tapped, polygons with N > 12 show a sharp transition 
with a clear density gap. On the other hand, triangles and squares present density 
fluctuations that are larger at large tapping intensities in contrast with all other shapes 
(including disks). 

A number of questions arise from this study that can lead future research. Some 
of these questions are: 

(i) What is the true nature of the transition for polygons with a large number of 
vertices? Can this transition be effectively found in infinite rough disks? Can the 
low density coexisting phase be related with the so called random close packing 

state pa manning. 

(ii) Given that fluctuations have a different trend, is the granular (configurational) 
temperature in the case of triangles and squares radically different from that of 
other shape packings? 

(iii) Given that for pentagons the fluctuations are equivalent for states at each side of 
4>min obtained with different T, which suggest that the states are equivalent, is the 
force moment tensor equivalent? 

Acknowledgments 

We thank Ana Maria Vidales and Irene Ippolito for valuable discussions. This work has 
been supported by CONICET and ANPCyT (Argentina). 

References 

[1] E. R. Nowak, J. B. Knight, M. Povinelli, H. M. Jaeger, and S. R. Nagel, Powder Technol. 94,79 
(1997). 

[2] L. A. Pugnaloni, M. Mizrahi, M. C. Carlevaro, F. Vericat, Phys. Rev. E 78, 051305 (2008). 



Steady state of tapped granular polygons. 



10 



L. A. Pugnaloni, D. Maza, I. Sanchez, P. A. Gago, J. Damas, I. Zuriguel. larXiv: 100 2.3264 (2010) 
A. D. Rosato, O. Dybenko, D. J. Horntrop, V. Ratnaswamy, L. Kondic, Phys. Rev. E 81, 061301 
(2010). 

I. C. Rankenburg, R. J. Zieve, Phys. Rev. E 63, 061303 (2001). 

F. X. Villarruel, B. E. Lauderdale, D. M. Mueth, H. M. Jaeger, Phys. Rev. E 61, 6914 (2000). 

G. Lumay, N. Vandewalle, Phys. Rev. E 74, 021301 (2006). 

M. Ramaioli, L. Pournin, Th. M. Liebling, Phys. Rev. E 76, 021304 (2007). 
R. Cruz Hidalgo, I. Zuriguel, D. Maza, I. Pagonabarraga, J. Stat. Mech. P06025 (2010). 
Y. Limon Duparcmeur, J. P. Troadec, A. Gervois, J. Phys. I (France) 7, 1181 (1997). 
M. Ammi, D. Bideau, J. P. Troadec, J. Phys. D 20, 424 (1987). 

C. Radin, J. Stat. Phys. 131, 567-573 (2008). 

D. Aristoff, C. Radin. larXiv:0909. 26081 f2009). 
Y. Jin, H. A. Makse. larXTv:1001. 52871 (2010). 

A. M. Vidales, L. A. Pugnaloni and I. Ippolito, Phys. Rev. E 77, 051305 (2008) 
A. M. Vidales, L. A. Pugnaloni and I. Ippolito, Gran. Matter 11, 53 (2009) 
J. A. Dijksman, M. van Hecke, Eur. Phys. Lett. 88, 44001 (2009). 
Box2D Physics Engine, www.box2d.org 

E. Catto, Rerative dynamics with temporal coherence (2005), available at 



http://box2d.googlecode.com/files/GDC2005_ErinCatto.zip (retrived on October 2010). 



L. A. Pugnaloni, M. G. Valluzzi, L. G. Valluzzi, Phys. Rev. E 73, 051302 (2006). 
R. Arevalo, D. Maza and L. A. Pugnaloni, Phys. Rev. E 74, 021303 (2006). 



M. D. Shattuck, arXiv:cond-mat/0610839 (2006). 

S. F. Edwards, R. B. S. Oakeshott, Physica A 157, 1080(1989). 

S. McNamara, P. Richard, S. Kiesgen de Richter, G. Le Car, R. Delannay Phys. Rev. E 80, 031301 
(2009). 

M. Schrdter, D. I. Goldman, H. L. Swinney, Phys. Rev. E. 71, 030301 (R) (2005). 
R. Blumenfeld, S. F. Edwards, J. Phys. Chem. B 113, 3981 (2009). 



